Structure formation and the origin of dark energy 
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Cosmological constant a.k.a. dark energy problem is considered to be one major challenge in 
modern cosmology. Here we present a model where large scale structure formation causes spatially- 
flat FRW universe to fragment into numerous 'FRW islands' surrounded by vacuum. We show that 
this mechanism can explain the origin of dark energy as well as the late time cosmic acceleration. 
This explanation of dark energy does not require any exotic matter source nor an extremely fine- 
tuned cosmological constant. This explanation is given within classical general relativity and relies 
on the fact that our universe has been undergoing structure formation since its recent past. 

PACS numbers: 95.36.+x,98.80.-k,98.80.Jk 

Several recent experimental observations [l| seem to 
strongly suggest that we live in the universe whose en- 
ergy budget is dominated by contribution from a mys- 
terious source which is otherwise invisible or missing in 
direct observations. This energy component is generally 
referred as dark energy. Furthermore, the dark energy 
component appears to have negative pressure. Result 
from supernova observations 0] that universe is under- 
going a recent acceleration seems to confirm such peculiar 
behavior of dark energy. 

Many attempts have been made in literature to under- 
stand the origin of dark energy (see [|| for some reviews). 
Arguably the most economical one is to introduce a non- 
zero cosmological constant in Einstein equation. How- 
ever, there are several conceptual difficulties associated 
with it. Firstly, experimentally required value of cosmo- 
logical constant turns out to be ~ 10 -123 in Planck unit. 
Such an extremely fine-tuned value of cosmological con- 
stant appears to defy any hope of possible explanation 
from a fundamental theory. Second conceptual problem 
is the so-called cosmic coincidence problem: why does en- 
ergy density due to cosmological constant which remains 
unchanged, become comparable to changing matter en- 
ergy density only at current epoch? 

There have been other attempts where it is anticipated 
that late time acceleration could be due to a dynamical 
scalar field that evolves in a suitably engineered potential 
0| , due to super-horizon perturbations || , from averag- 
ing of LTB models @ or due to back-reaction effects of 
inhomogeneous structures in the universe @, S, @|- In 
particular, using Buchert equations one can get accel- 
eration by incorporating inhomogeneous back-reaction. 
However, same back-reaction term contributes negatively 
to energy density thus fails to explain missing energy 
unless one also introduces strong negative spatial curva- 
ture. In brief, none of the current attempts can explain 
the origin of dark energy satisfactorily. 

According to standard model of cosmology, continued 
expansion of the universe causes radiation to eventu- 
ally decouple from matter. Subsequently, largely homo- 
geneous matter distribution with small inhomogeneity 



FIG. 1: A simple illustration of late time universe where 
structure formation causes FRW universe to fragment into 
numerous 'FRW islands'. Each FRW island (shaded region) 
is surrounded by vacuum metric. 



starts collapsing to give rise current large scale struc- 
tures. Initial phase of structure formation can be de- 
scribed by linear perturbation theory around homoge- 
neous background. However, in later phase when struc- 
ture formation enters non-linear regime such descriptions 
are insufficient. In this era, matter distribution consists 
of numerous dense regions surrounded by relatively rarer 
regions. For simplicity here we make sharp-boundary ap- 
proximation for the matter distribution during later pe- 
riod of structure formation such that matter is contained 
within spherical, homogeneous regions surrounded by 
empty space. Boundary of such spherical regions, how- 
ever, needs to be continuously refined as on-going struc- 
ture formation continues to cause dense regions to be- 
come denser. In other words, in this model structure for- 
mation along with sharp-boundary approximation causes 
FRW universe to fragment into numerous 'FRW islands' 
with shrinking boundary that are surrounded by vacuum 
(see FigfT|). We will show that this mechanism can ex- 
plain the origin of dark energy as well as the late time 
cosmic acceleration. 

In the standard model of cosmology, spacetime describ- 
ing our universe is assumed to be foliated by homoge- 
neous and isotropic spatial hyper-surfaces parametrized 
by a global time. We imagine such an observer who treats 
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FIG. 2: A spherical patch of universe of coordinate diameter 
L which was homogeneous at the beginning of structure for- 
mation. In a later period the matter distribution, by means of 
sharp-boundary approximation, is contained within a spheri- 
cal region (shaded) of coordinate diameter U. 

the spatial hyper-surfaces as homogeneous and isotropic 
during entire evolution of universe and accordingly mea- 
sures distance using an average spatially-flat Friedman- 
Robertson- Walker (FRW) metric 

ds 2 = -dt 2 +a 2 (t)dx 2 , (1) 

where a(t) is the scale factor. We refer this observer as 
observer A. We consider another observer, say observer 
B, who uses same time parametrization of the spatial 
hyper-surfaces as observer A, but is careful to consider 
the effects of structure formation. In particular, to mea- 
sure distance during structure formation observer B uses 
flat FRW metric g % — diag{— 1, af , af , af) inside the 
spherical regions containing homogeneous matter distri- 
bution otherwise uses vacuum metric g° . To illustrate 
this, let's consider a spherical patch of universe of coor- 
dinate (also co-moving for observer A) diameter L which 
is homogeneous at the beginning of structure formation 
(see FigUJ. The patch then begins to undergo structure 
formation such that at a given time collapsing matter dis- 
tribution, by means of sharp-boundary approximation, is 
contained within a spherical region of coordinate diam- 
eter li(t) with < U < L. We assume that the vacuum 
metric g° v describing the spherically symmetric empty 
space created due to structure formation, is static. (This 
would be the case if we ignore the presence of neigh- 
bouring patches, as then Birkhoff's theorem would imply 
spherically symmetric vacuum metric is static.) The vac- 
uum metric g°„ however is inhomogeneous as can be seen 
from continuity of metric at the boundary of inside re- 
gion. 

To have a precise notion we define average FRW met- 
ric ((T|) such that the proper distance between the points 
Pi and P2 (see FigJ5]) as measured by observer A is equal 



to the proper distance measured by observer B. With- 
out loss of generality we assume that the points lie on 
a^-axis. Given elementary proper distance in a spatial 
hyper-surface is | y/ gjkdx^dx k where j, k represent spa- 
tial coordinates, the definition leads to 



aL = ai l t + 2 Jgf 1 (x 1 ,x 2 ,x 3 )dx 1 . (2) 

hill 

Using the definition, we can also compute relation be- 
tween expansion rates of the points as measured by both 
observers. In particular, 

aL = — / ^/ gn{x 1 ,x 2 ,x 3 )dx l 

dt J-L/2 

= aJi + aih - li^Jg^ik/ 2,0, 0) = dik , (3) 

where over-dot denotes derivative w.r.t. to time t. In 
last line, we have used continuity of metric solution i.e. 
g° 1 (l i /2,Q,0) — g\ x = af. We may note that to have 
an expanding average metric, the metric of the inside 
region must also be expanding. However, proper volume 
of the inside region can still decrease as its coordinate 
diameter shrinks. We now define average energy density 
p for observer A, by requiring that at any given time total 
energy contained within the patch of coordinate diameter 
L is same as measured by observer B. In particular, if 
observer B measures energy density of inside region to 
be pi and of outside region to be zero then 

o=(S) 3 <" < 4 > 

As defined, n 3 is the fraction of total proper volume oc- 
cupied by matter distribution. Physically, parameter n is 
a measure of amount of structure formation and satisfies 
< n < 1. In particular, n = 1 implies there are no 
underlying structures. This can be seen from the equa- 
tion |[2}. As we will see, equations ([3]) and (j4]) form the 
backbone of arguments presented here. 

Let's imagine that both observers want to confront 
their respective Einstein equations with experimental 
data. Observer A performs separate observations to mea- 
sure expansion rate as well as energy density. However, it 
turns out that to make a right balance, observer A needs 
to postulate an extra invisible component in Friedmann 
equation for the average metric i.e. 

3 (£j =8irG{p + p DE ) , (5) 

where G is Newton's constant and pde denotes dark en- 
ergy component. On the other hand observer B uses 
standard Friedmann equation for the region containing 
homogeneous matter distribution and uses vacuum Ein- 
stein equation for the remaining region. In particular, 
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Friedmann equation for observer B is 
:-! | — | = 8nG Pi . 



(6) 



Equipped with the details of underlying structures i.e. 
using equations ([3]), and Friedmann equation |6]), ob- 
server B can derive Friedmann equation for average FRW 
metric ([T]), given by 



(7) 



One may note that right hand side of equation (J7]) has 
an extra energy density component apart from average 
energy density p. Thus, comparing equations ([5]) and (J7J), 
observer B can derive the expression of dark energy that 
observer A should perceive 
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(8) 



In the situation when n = 1, dark energy component dis- 
appears. In other words, observer A wouldn't have per- 
ceived any dark energy component if there were no under- 
lying structures in universe. Since existence of underlying 
structures requires parameter values to be < n < 1 then 
dark energy component is necessarily positive. This is in 
contrast with back-reaction models such as 0, Q where 
back-reaction term contributes negatively to energy den- 
sity. As evident, dark energy component is comparable 
to the magnitude of average energy density p. So this 
model can naturally explain cosmic coincidence prob- 
lem. Another crucial property of dark energy component 
([8]) is that it may appear as constant even though it is 
naively proportional to decreasing average energy density 
p. With the beginning of structure formation, the value 
of parameter n starts decreasing from unity. So during 
structure formation proportionality factor (1/n — 1) in- 
creases. This implies that for suitable rate of structure 
formation dark energy component (jHJ) may appear as con- 
stant. 

For observer A total energy contained within the patch 
is given by E = V p where proper volume of the patch 
V = (ira 3 L 3 /6). Using the definition of pressure P — 
— (dE/dV), one can derive the conservation equation for 
energy density p as 

,=-(;)(Mf) = -^u+.>, e, 

where H := (a/a) is Hubble parameter and us := P/p is 
the corresponding equation of state. Analogously, we can 
define equation of state u>de for dark energy component 
such that pde = —3H pde(1 +^de)- Using equations 
© and ©, we can compute equation of state for the 
dark energy 



where (h/n) =: —r n H. Parameter r n is a measure of 
rate of structure formation. Dark energy expression (jSJ) 
along with its equation of state (|10[) can mimic a cosmo- 
logical constant at current epoch for suitable values of 
structure formation parameters n and r n . In particular, 
the values of structure formation parameters such that 
r n = 3(1 — n)(l + us) will lead to uide = —1 which is 
the equation of state for a cosmological constant. How- 
ever, a distinguishing feature between them is that while 
equation of state for a cosmological constant remains un- 
changed, the equation of state (|10[) varies with time. 

To explicitly show that structure formation can lead to 
an accelerating phase, it is convenient to compute Ray- 
chaudhuri equation. Taking time-derivative of equation 
([7]) and then using conservation equation ([9]) , we can de- 
rive Raychaudhuri equation as 



3|-|= 4ttG 
a 



r^ _ (1 + 3us) 
n 



n 



P 



(11) 



From equation (fTT|) it can be seen that the values of struc- 
ture formation parameter such that r n > (1 + 3a;) will 
lead to an accelerating phase for observer A. 

The modifications to average dynamics of the given 
patch due to underlying structures, depend on the values 
of parameters n and r n . However, the modifications do 
not depend explicitly on coordinate diameters L and 
Thus, if one considers a different patch but with same 
values of parameters n and r n , then one will get same 
Friedmann and Raychaudhuri equations. Given one can 
pack R 3 space very closely using 3-spheres of arbitrary 
diameters hence the modified Friedman equation ([7]) and 
modified Raychaudhuri equation (fTTj) can be considered 
as good approximation of equations that describe average 
dynamics of the model universe with underlying struc- 
tures. Observed contribution from the dark energy com- 
ponent at current epoch is about 70% of the critical en- 
ergy density. The equation {8} with the parameter value 
n = 0.3 can lead to such an observed amount of dark 
energy. If we consider the average matter to be pressure- 
less i.e. us = then lode = — 1 requires parameter value 
r„ = 2.1. 

To characterize the matter distribution of inside re- 
gion as seen by observer B, we need to compute rela- 
tion between average equation of state us and equation 
of state for inside region uji := Pi/ pi where pressure 
Pi = —(dEi/dVi). Ei and Vi are total energy and proper 
volume of the inside region respectively. As earlier we 
can derive the conservation equation for inside region, 
given by 



usde = -1 



(l + o,)- 



3(1 -n)_ 



(10) 




(12) 
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Coordinate diameter of inside region is time-dependent 
and it is reflected in conservation equation (|12p with its 
explicit dependence on (li/li). Using the relation between 
energy density ([3]), their conservation equations © and 
(|f 2p . one can compute relation between the equation of 
states 



w = (1 - r n )uji 



(13) 



For gravitational collapse to occur with ordinary matter, 
the corresponding matter distribution must be pressure- 
less. So for observer B, matter distribution of inside re- 
gion should be pressure- less i.e. LUi = 0. On the other 
hand observer A finds average matter also to be pressure- 
less i.e. U — 0. The equation (fi~3|) ensures that physical 
requirement of observer B and observed fact for observer 
A can be consistently met. We should mention here that 
to derive equation (fTTj) . one can also use Raychaudhuri 
equation for observer B. However, in doing so one should 
be careful to include the additional pressure component 
coming from the shrinking boundary. 

Experimental observations seem to also imply that to- 
tal energy of the universe has another dark component, 
the so-called dark matter which is pressure-less. From 
equation (|10p . one may note that if u> = and r n = 
i.e. if (—li/li) = (1/n — 1)H then the equation of state 
mimics a pressure-less energy component. It is conceiv- 
able that some 'FRW islands' may have different values of 
structure formation parameters leading to such behavior. 
However, whether such scenario can explain phenomena 
ascribed to the presence of dark matter such as galaxy 
rotation curves, remains to be explored. 

To summarize, we have argued that the origin of dark 
energy can be understood as a consequence of large scale 
structure formation. This explanation of dark energy 
does not require any exotic matter source nor a fine- 
tuned cosmological constant. However, presented model 
in its current form has several deficiencies. Firstly, we 
assume that structure formation leads to creation of void 
around each FRW island. However, we know cosmic mi- 
crowave background (CMB) photons fills up entire uni- 
verse. Thus, even though CMB contribution to average 
energy density is negligible during structure formation 
but for an accurate description one should consider their 
presence. In this model net effects of structure formation 
on dynamics of average metric can be summarized by in- 
troducing just two characteristic parameters n and r n . 
However, the model itself does not shed any light on the 
values of the parameters n and r„. We may recall that 
the model is based on sharp-boundary approximation of 
the matter distribution which is under-going structure 
formation. Thus, to compute relation between the val- 
ues of the parameters one needs to perform a detailed 
simulation of structure formation with a matter distri- 
bution which should be then successively approximated 



by sharp-boundary approximation. Finally, we have not 
addressed the issue: why is cosmological constant zero in 



our universe.' 
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